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The number of modes that can propagate without radiation loss in 
oversized waveguides is sharply reduced if the waveguide is coupled to a 
structure supporting radiation modes, the loss mechanism being analogous 
to Cerenkov radiation. The coupling formula derived in Part I 1 is used 
to evaluate the loss for a specific configuration: a reactive surface (e.g., 
a thin dielectric slab) acting as a waveguide, coupled to a semi-infinite 
dielectric acting as a mode si?ik. The method consists in first assuming 
that the substrate is finite in size and lossy and adding the losses associated 
with each substrate mode. The substrate dimensions are subsequently 
made infinite and the dissipation loss is made to vanish. The expression 
obtained for the radiation loss coincides with an expression obtained by 
solving the boundary value problem. The method is then applied to the 
problem of mode selection for dielectric rods coupled to dielectric slabs, 
which is of particular importance for optical communications and inte- 
grated optics. A 2-dB/m radiation loss is calcxdated for the first higher 
order mode whei\ the rod radius is 10 pm, X = 1 fim, n = 1.41, and the 
rod-to-slab spacing is 0.15 pm. 

I. INTRODUCTION 

An expression for the coupling between lossy single-mode open wave- 
guides was derived in Part I. 1 We now investigate the coupling of a 
waveguide with finite cross section with a waveguide with infinite 
cross section (called a substrate), the latter supporting radiation modes. 
Radiation losses are suffered whenever the propagation constant h of 
the guided mode is smaller than the highest propagation constant h. 
of the radiation modes carried by the substrate. Radiation then takes 
place at the Cerenkov angle = cos -1 (h/h,). By properly choosing 
the dimensions and permittivities of the waveguide and those of the 
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substrate, it is possible to reduce the number of modes that can 
propagate without attenuation (in the absence of dissipation and 
scattering losses). This arrangement is of great practical importance 
because optical fibers are usually highly overmoded to facilitate fabrica- 
tion and splicing. 2 (For a coherent source, it is important to reduce the 
number of modes because different modes usually have different 
group velocities. If a short optical pulse is sent through the fiber, mode 
conversion takes place because of the imperfections of the fiber; this 
causes the pulse to spread in time.) The mode selection mechanism 
just described is also of practical importance in the microwave range 
for oversized waveguides such as oversized microstrips on dielectric 
substrates and oversized dielectric strips.' Multimoding in traveling 
wave tubes can also be avoided with the help of mode sinks. 

We investigate the loss mechanism for two specific configurations. 
First, a reactive surface acting as a waveguide coupled to a semi- 
infinite dielectric acting as a mode sink. We show that, by adding the 
losses associated with each substrate mode, an expression for the total 
loss is obtained that coincides with an expression obtained by solving 
the boundary value problem. Then the method is applied to the prob- 
lem of a dielectric rod coupled to a dielectric slab. 2 The case of dielectric 
rods coupled to dielectric cylinders supporting whispering gallery 
modes and acting as mode sinks 3 will be discussed in another paper. 

II. RADIATION LOSSES IN SUBSTRATES— GENERAL FORMULA 

To evaluate the radiation losses, let us first assume that the trans- 
verse dimensions of the substrate are finite, and let h lt = h tr + ih ai 
be the propagation constant of a trapped mode in the substrate, with 
h tr real and h ti real positive (the subscript s stands for "substrate"). 1 
If h B denotes the propagation constant of a trapped mode of the wave- 
guide in the absence of the substrate, the propagation constant h of 
the coupled wave is, from eq. (6a) in Part I, 

h - h„ + h(h„ - K) - [\{h» - ho) 2 + C 2 ?, (1) 

where C 2 = c a c h /P a Pb denotes the coupling coefficient defined in Part 
I. The minus sign before the square root has been selected because it 



*In the microwave range, there are no compelling reasons for using dielectric 
waveguides that are large compared with the wavelength in all dimensions, but we 
may want to use strips (either metallic or dielectric) whose widths exceed one 
wavelength for improved accuracy. . , 

t The dependence of the field on time (t) and on the axial coordinate (z) is denoted 
exp [i(hz — M03' This term is henceforth omitted. 
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corresponds to the mode whose field is concentrated in the waveguide 
cross section rather than in the substrate (that is, we require h = h 
when C 2 = 0). 

Let us now assume that h is real (lossless waveguide) and that 

h.i » C. (2) 

Using this condition, eq. (2), we can expand the r.h.s. of eq. (1) in 
power series of C 2 and keep only the first two terms in the expansion. 
The loss is given by the imaginary part hi of h. Because the imaginary 
part of C 2 can be neglected in the case that we consider, we have 

hi « C*h.it(h„ - KY + hlJ-K (3) 

The total loss £ experienced by the waveguide is now obtained by 
summing over the various modes of the substrate : 

£ = Z CVi.Uh.ra - h o y + /&]->, (4) 

where the subscript a refers to the substrate modes. We have assumed, 
for simplicity, that h ti does not depend on a. It is shown in the next 
section for a simple configuration that in the limit of dense substrate 
modes eq. (4) is in agreement with an exact result, obtained from a 
boundary value method. 

If we let the cross-section area S of the substrate tend to infinity, 
the substrate modes become denser and denser, and the summation 
in eq. (4) can be replaced by an integral 

£ = lim £ C 2 a h.l(h. ra - h„y + hIJr 1 

S-»°o a 

= J e(h. r )h.l(k. r - h y + h 2 H -]- l dh. r , (5) 

where we have defined a coupling density e by 
e(h. r )dh. r = lim £ C% 

S-oo a 

the range of a being defined by the condition 

h, r < h, Ta < h, r + dh, T . (6) 

This density exists because, as S — » °o , the coupling coefficient C 2 
decreases at least as fast as S~ l , the power in the substrate being 
proportional to S if the power density is kept a constant. 

We can now let h.i tend to zero, the condition eq. (2) being pre- 
served. The second factor in the integrand of eq. (5) is sharply peaked 
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at h ar = K and behaves as a symbolic 5-function. Thus, in the limit 
h,i —* we have 

£ = ire (ho). (7) 

It should be noted that the subscript a in eqs. (4) to (6) stands for 
three subscripts m, n, and s, where m refers to modes in the x direction, 
n refers to modes in the y direction (we assume for simplicity that the 
substrate modes are separable in Cartesian coordinates), and s refers 
to the state of polarization (e.g., H or E modes). 

III. COUPLING TO A SEMI-INFINITE SUBSTRATE 

Consider a reactive surface coupled to a semi-infinite dielectric 
(Fig. 1). We consider only H modes and assume that the field is 




Fig. 1 — (a) Reactive surface, with normalized susceptance a, coupled to a semi- 
iufinite dielectric with permittivity « = n 2 e„. For H modes, the structure is assumed 
terminated in the y direction by electric walls. Radiation takes place at the Cerenkov 
angle = cos -1 t(k* + a»)*/*»]i k = 2"" A- (°) Variation of the field as a function 
of x. 
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independent of the y coordinate. Except for the changes x <=* y and 
i — > — i, we use Shevchenko's notation. 4 

For waves propagating along the z axis, the electric field has only a 
y component that we denote E. In a region with constant e, E obeys 
the wave equation 

d 2 E/dx 2 + (w 2 eM„ - h 2 )E = 0, 

H x = -hi^o^E, 

Hz = (wito^dE/dx. (8) 

If e has a finite discontinuity, E and dE/dx remain continuous. 
The general solution of eq. (8) for e = e and e are, respectively, 

E = A+e'** 4- A~e-** m t (e ) (9a) 

E, = A+e*' + Are-<°*, (e) (9b) 

X 2 = co 2 € oMo - h 2 , (10a) 

g 2 = co 2 en -h 2 = u 2 + X 2 , 
W z s w 2 (e - c )m - 

The loss can be evaluated by solving the boundary value problem . 
At the reactive surface (x = — D), we have the condition (see Ref. 4) 

dE/dx + aE = 0, x = -D, (11) 

where a is a positive real number proportional to the susceptance of the 
surface. * We assume that, in the dielectric, the wave propagates away 
from the structure, that is, 

E, = Af&*. (12) 

Note that h is expected to have a small positive imaginary part ex- 
pressing the radiation loss in the dielectric. Assuming that e is real, 
that is, that the dielectric is free of dissipation losses, eq. (10b) shows 
that g has a small negative imaginary part. Thus, the wave amplitude 
grows exponentially as the distance to the structure increases. This 
solution of Maxwell's equations is called a "leaky wave." 4 It is not 
difficult to show that the curves of constant irradiance in the dielectric 
are straight lines making with the z axis an angle 9 = cos -1 (h /kn) 
(Cerenkov angle). 



* A thin dielectric slab with permittivity t and thickness d, supported by a 
magnetic wall, is equivalent to a reactive surface with normalized susceptance 
a = « s (e — t )nod. An equivalent configuration, obtained by symmetry with respect 
to the magnetic wall, is a thin slab of width 2d with dielectrics symmetrically located 
on both sides. Note that a has the dimension of a propagation constant. 
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From eq. (12), the boundary condition at x = is 

dE/dx - igE = 0, x = 0. (13) 

From eqs. (9a), (11), and (13), we obtain the equation defining X, or h, 

(X - ia)(X + fir) = (X + ia)(X - g) exp (2ftD). (14) 

If we let aD tend to infinity, the reactive surface is uncoupled to the 
dielectric and eq. (14) reduces to X = X = ia; that is, 

X 2 = X 2 m u>*e Ho - h\ = -a 2 , (15a) 

g * = g * = fa ,2 (e _ eo ) Mo + xj. (15b) 

Equation (15a) defines the propagation constant h of the uncoupled 
reactive surface. 

Let us now consider 

exp (2iX D) m 5 (16) 

as a small parameter and set 

X = X + Xx5 + • • •, 

g = go + gi& + • • • > 

in eqs. (14) and (10b). Collecting terms of first order in 5 we get 
X! = 2ia(ia - g,)/(ia + g„). (18) 

From eqs. (10a) and (17) we have, to first order, 

Im(h) = -(a8/h )Re(X 1 ). (19) 

Thus the loss £ ■ lm(h) is 

£ = 4a 3 w" 2 /i - Vo exp ( - 2aZ)) , (20a) 

or, explicitly, in terms of A;, n, D, and a, 

£ = 4a 3 [fc 2 (n 2 - l)]" 1 ^ 2 + <* 2 ) -i 

X [fc 2 (« 2 - 1) - « 2 ] § exp (-2oD). (20b) 

If the micron is used as the unit of length, the loss in dB/km is ob- 
tained by multiplying the r.h.s. of eq. (20b) by 8.7 X 10 9 . 

This expression for the loss, applicable to small couplings, can be 
obtained alternatively from the equality 

h- h = u /"(« - e )E+-E p dS //(E + X H p - E p X H + )-dS, (21) 
where (E, H) and h denote the field and propagation constant of the 
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wave guided by the reactive surface in the absence of the dielectric 
and (E+ H+) denotes the field adjoint to (E, H) (see Part I). (E p , H„) 
and h denote the field and propagation constant in the presence of the 
dielectric. The integral in the numerator extends to the dielectric cross 
section, and the integral in the denominator extends to the whole 
cross section. Equation (21) is exact and is readily obtained from Max- 
well's equations.* The field (E p , H p ), unfortunately, is not known. It 
may differ considerably from the unperturbed field (E, H) when the 
dielectric supports modes almost synchronous with the waveguide 
mode. This is why this expression, eq. (21), is, in general, not practical 
to evaluate the coupling between waveguides, or waveguides and mode 
sinks. The configuration presently considered, however, is sufficiently 
simple to be handled on the basis of eq. (21). 

For our case, eq. (21) becomes, with the approximation h £b h , 

h- h tt - (u>W2/i„) r (e - e )EE p dx/r EHx. (22) 

The unperturbed field, normalized to unity at x = — D, is 

E = exp (ixx) exp (ixD). (23) 

The perturbed field is obtained by assuming as before an exp (igx) 
dependence in the dielectric, matching E and dE/dx at the vacuum- 
dielectric interface (x = 0), and stating that E p £Hl at x = — D. 
We obtain 

E p = 2(1 + 0/X)- 1 exp (ixD) exp (igx), x ^ 0. (24) 

Substituting in eq. (22) and integrating, a result identical to eq. (20) 
is obtained. 

Let us now apply to the same problem the method explained in 
Section II of this paper, which consists in adding the losses associated 
with each mode of the substrate. The coupling coefficient between two 
H modes, with fields E and E„, was given in Part I. With our present 
notation we have 

C 2 = aViH E 2 / Je^x^E 2 ,/ fE]dx\ (25) 

where the integrals are over the whole cross section, and E, E, are 
defined at some point located between the two waveguides. 



* The contribution at infinity is assumed to vanish. Thus, it is implicitly assumed 
that the rate of decay of the unperturbed field exceeds the rate of growth of the 
perturbed field. This condition is always satisfied for small couplings. 
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The field E of the reactive surface alone is, as we have seen, 

E = exp (-ax). (26) 

Thus, at x = 0, 

E 2 /r EHx = 2a exp (-2aD). (27) 

Let us consider next the dielectric alone and first assume that its 
thickness L x is finite. By matching E and dE/dx at x = and x = L x , 
we obtain the field at the vacuum dielectric interface, and 

&,/ /" +e ° E\dx « E\ /J*' E\dx = 2glir>L x - 1 . (28) 

Substituting eqs. (27) and (28) in eq. (25), we obtain the coupling 
coefficient 

C 2 = 4a 3 glu- 2 K 2 exp {-2aD)L^. (29) 

Let us now evaluate the number of modes (N dh) in the dielectric 
whose propagation constants lie between h and h + dh. Because we 
are far from cut-off, the boundary condition is almost the same as for 
a metallic waveguide, E = 0. Thus, the condition on g is 

g m = mr/L x , m = 1, 2, ••-. (30) 

Using the relation 

gl = coVo - h\ (31) 

the mode number density is, from eq. (30), 

N = hg-iLJ*. (32) 

The radiation loss is obtained from eqs. (29), (32), and (7), and h = h , 

= Qo, 

£ = V C 2 N = ±c?u- 2 K l g exp (-2aD). (33) 

This result coincides with the result eq. (20) obtained by taking the 
limit of large D in the exact solution. The variation of the loss expressed 
in dB/km is given in Fig. 2 as a function of the normalized susceptance 
a of the surface, for X = 1 /mi, e/e = 2, and D « 1.5, 1.75, and 2 /tm. 
For comparison, when the dielectric permittivity has the form 
e = €o -f- { €i (the dielectric is perhaps a lossy foam) and the spacing 
D is chosen as large as consistent with a loss of 10 dB/km at a = 6.28, 
the loss experienced is shown on the same figure as a dotted line. The 
comparison clearly shows the advantage of mode sinking over dis- 
sipation for mode selection. 
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Fig. 2 — Radiation loss in dB/km as a function of the normalized surface suscep- 
tance a of the waveguide for a wavelength X = 1 nm, n 2 = 2, and D = 1.5, 1.75, 
and 2jim. The dotted line is applicable to a dissipative dielectric. 

IV. COUPLING TO PLANAR SUBSTRATES 

Let us now consider a waveguide with propagation constant hi 
coupled to a substrate that extends to infinity in the y direction, but 
has a finite thickness in the x direction. This substrate is perhaps a 
reactive plane (e.g., a corrugated conductor) or a dielectric slab, as 
illustrated in Fig. 3. In any case, homogeneity of the substrate in the 
y, z plane is assumed. 

Because of the assumed homogeneity of the substrate, plane wave 
solutions 

E.(x, y, z) = E a (x) exp (ih, v y + ih tt z), (34) 

where 

h„ = f(h, u , o)), (35) 

exist at some angular frequency u> (ui is now considered a fixed param- 
eter and is omitted). 



TRANSVERSE COUPLING IN FIBER OPTICS 



683 




Fig. 3 — Dielectric rod coupled to a dielectric slab. The rod field E is shown for the 
spurious Hoi mode, and the slab mode is Hi, (p is a continuous index in the limit 
L v — » oo ). Coupling takes place at 4> » 0. 

In the discussion that follows, we consider only waveguide and sub- 
strate modes that are even in y. Assuming that / is even in h, y and 
that the slab is terminated by electric walls, even modes satisfy the 
relation 

h. y L y = 2mr, » = 0, 1, 2-«, (36) 

where L„ denotes the width of the substrate. L v will be later assumed 

to tend to infinity. The density N of even modes is from eqs. (35) 

and (36) 

N = (df/dhJ-^Ly^ir). (37) 

If the substrate is isotropic, with wave vector h„ eq. (35) is 

h„ ^ f(h. v ) = (ti - /&)*, (38) 

and the mode density is, from eq. (37), 

N = {hJh. v )L v /2T. (39) 
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The loss is then obtained from eqs. (7) and (39). 

£ = HVf-W]^ (40) 

the coupling coefficient C 2 being evaluated from eq. (6) in Part I. 

It should be noted that, when the propagation constant of the wave- 
guide mode (hi) is just equal to the propagation constant (h.) of the 
2-dimensional substrate, h,„ is equal to zero and the loss, according to 
eq. (40), is infinite if C 2 L„ remains finite. (This was not the case for 
the 3-dimensional mode sinks considered in Section III because, as 
L x — > «> , the field at the surface of the dielectric tends sufficiently 
rapidly to zero to make C 2 L X vanish in the limit.) This infinity at 
hi = h a would be removed if some finite dissipation loss in the substrate 
were present. Even in the absence of dissipation losses, the radiation 
loss remains finite at hi = h 3 , because the perturbation method on 
which eq. (40) is based is no longer applicable. The peak in the loss 
curve predicted by eq. (40) (analogous to a sound barrier) is pro- 
nounced only for small couplings. 

Our general result, eq. (40), is now applied to a dielectric rod coupled 
to a dielectric slab. The thickness and permittivity of the slab can 
always be chosen in such a way that only the fundamental mode of the 
rod propagates without radiation loss. The calculation of the loss of 
higher-order modes is carried out for the case where the rod diameter 
and the slab thickness are very large compared with the wavelength; 
that is, when the rod is highly multimoded in the absence of coupling. 

Approximate expressions for the modes and propagation constant 
in the slab and the rod are given in the next subsections. 

4.1 Modes of the slab 

Let us consider first the modes in the dielectric slab. If the thickness 
2d of the slab is very large (more precisely, if co 2 (e — e )f*od 2 » 1), the 
propagation constant of the fundamental Hi mode is approximately 
given by the condition that the field E vanishes at the boundary 

E(x, z) tt E so cos (g a x) exp (ih a z). 

Thus, we have* 

g 2 s = o> 2 e/z - h 2 s = (w/2dy. (41) 



* A more accurate and general expression is (see Ref. 5) gd = m(ir/2){\ — V~ l ) 
for H modes and g,d = m(ir/2)(l — n~ 2 F -1 ) for E modes, where m = 1, 2- • • is the 
mode number and V = ud. These expressions show that the Hi mode that we are 
considering in this section is the fundamental mode ; that is, the mode that has the 
largest propagation constant. The difference AA in propagation constants is, for 
m = 1, equal to dr l (w/2kndY{\ - 1/n 2 )*. 
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(The axial coordinate is denoted z instead of z to avoid changing our 
notation when waves propagating at some angle to the z axis are con- 
sidered. The origin of the x axis is, in this subsection, at the center 
of the slab.) The axial (z') and transverse (x) components of the 
magnetic field are, within the slab, as we have seen before 

H x = -h,(uti )- l E, (42) 

H t . = {wno)~ l dE/dx, (43) 

and the power per unit width is approximately 

Ptt-f +d EH x dx = dh.(uiio)-*El. (44) 

The field at the boundary is in fact not exactly equal to zero. To 
obtain its value, we use the fact that the dependence of E on x in 
vacuum is exp ( — p a x), where p\ = h 2 , — w 2 e Ho, and the continuity of 
dE/dx. We obtain 

E(d) = (ir/2d)p^E, . (45) 

Now let the slab have a finite width L y with electric walls at 
y = ± Ly/2. The modes even in y can be described as a superposition 
of two infinite slab waves whose propagation constants are such that 

h ly = ± 2irn/L y> n = 0, 1, 2, • • -. (46) 

We have, by definition, 

k% + hi = h 2 s> (47) 

h s being given in eq. (41). 

The field has all its components different from zero with the exception 
of E x , which vanishes. The components E y and H z are obtained by 
adding the field of the two waves. We obtain 

E sv = 2hJh l cos Qi, y y) cos (irx/2d)E S0 , (48) 

H, z = -2(iwn )- 1 h, z hr 1 (ir/2d) cos (h. v y) sin (irx/2d)E, . (49) 

The energy flowing through the slab is obtained by multiplying P, 
given in eq. (44), by 2h, z h^ 1 L u 

P a = 2h tz {^ )-HL y El. (50) 

The y component of the field at the boundary (x = d) is obtained 
from eq. (45) or directly from H az = (iu)fi )~ l dE/dx: 

E ty (d) = -p s -Hun H sz . (51) 
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4.2 Modes of the rod 

Let us now turn our attention to the modes of the dielectric rod. 
We assume that the radius a of the rod is much larger than the wave- 
length (a» X). 

In the limit of large radii, the propagation constant of the funda- 
mental HEn mode is given (see the appendix) by the first root of 
Jo(g a), namely, 

g a = (wVo - A5)*o = 2.4- • •, a -> <». (52) 

The next higher order mode of the dielectric rod is the H Ql mode.* 
In the limit of large radii, the boundary condition at r = a is E$ = 0, 
as for a round metallic pipe. The propagation constant hi is therefore 
given by 

Ji(gia) = 0, (53) 

whose first root is 

gia = (co 2 e Mo - h\Ya = 3.8- • •, a -> *>. (54) 

Within our approximation, the field of the H i mode in the rod (r < a) 
has components 

E+ = Jiigir), 

H r = -h 1 (o )fio )- 1 J 1 (g 1 r), (55) 

H z = (iufi )~ l giJo(gir), 
and the energy flow is 



P = 



= - [* E+H r 2irrdr = *fti(6W»)~~Wo(0ia)- (56) 

To obtain the field E<, at the boundary (r = a), we use the fact that 
dE/dr is continuous and that the r dependence of E^ in vacuum is 
approximately 1 exp ( — pir) where p\ = h\ — co 2 e ju . We obtain 

E*(a) = pr l ianoH t . (57) 

4.3 Synchronization conditions 

For simplicity and because this is a case of practical significance, 
we assume that the rod and the slab have the same permittivity e. 



* The Em and HE21 modes have almost the same propagation constant as the Hoi 
mode for large rod radii. For small radiation losses, they can be considered indepen- 
dently of the Hoi mode (see appendix). 

t The exact dependence of E$ on r is Ko(pir), where Ko denotes the modified 
Bessel function of the second kind. For large arguments, Ko(x) » (2/irx)* exp (— x) 
and K' a (x) « - (2/tx)» exp (-x) ~-K (x). 
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The fundamental HEu mode of the rod is free of radiation loss if its 
propagation constant h„ given in eq. (52) is slightly larger than the 
propagation constant h, in the slab given in eq. (41). For simplicity, 
we set h, = h or, equivalently, g, = g ', that is, 

*/2d = 2.4/o, (58a) 

or 

d = 0.65a. (58b) 

Thus the ratio of the slab thickness to rod diameter is 0.65. (In practice, 
the slab has finite dissipation losses and a finite width. Furthermore, 
it is difficult to control accurately the thickness of the slab. For these 
reasons, it might be preferable to choose the value of h, midway 
between the propagation constants of the HEu and Hoi modes rather 
than equal to the propagation constant of the HEu mode. If the 
former condition were to hold, we would find that the slab thickness 
should be equal to half the rod diameter.) Figure 4 gives the propa- 
gation constants of the rod and the slab for n = 1.41 and a rod radius 
of 10 nva (X = 1 nm). 

Let us now consider one of the next higher order modes of the rod, 
the Hoi mode. This mode radiates into the substrate modes that have 
the same propagation constant along the z axis {h az — hi). Using eq. 
(54), we obtain 

« 2 € Mo - k% = (3.8/a) 2 . (59) 

Since 

h% + h% = h 2 s , (60) 

and h, has the value h given in eq. (52), we have 

h% = (3.8/a) 2 - (2.4/a) 2 , (61) 

or 

h av = 3.0/o. (62) 

In the next subsection we evaluate the coupling coefficient between 
the Hoi mode of the rod and the substrate mode defined by eq. (62). 

4.4 Coupling coefficient 

The contour of integration for the evaluation of the coupling 
coefficient being arbitrary, it is convenient to choose this contour as 
the rod boundary, r = a. Along that contour, the Hoi mode field is a 
constant. 
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Fig. 4 — Propagation constants (h) of the trapped modes of the rod and maximum 
value (h,) of the propagation constants of the radiation modes in the slab. It is 
assumed that n = 1.41, X = 1 /xm, and a = 10 /xm. The modes circled are those 
whose coupling is discussed in this paper. 



Let 4> denote the angle from the x axis shown in Fig. 3 and D the 
spacing between the rod and the slab. We have 



x = —D — a(l — cos <t>), 
y = a sin 0. 



(63) 



Because a » X, the coupling takes place near the point of closest 
approach of the rod to the slab ; that is, K, 0. We can therefore write 



xtt -D - a<f> 2 /2, 
y ^ a<j). 



(64) 



The y dependence of the field slab is cos (h, u y) = cos (h. y a4>). How- 
ever, since, according to eq. (62), h, y is of the order of ar\ the argu- 
ment of the cosine function is small compared with unity in the range 
where the coupling is significant. Thus, we can neglect the dependence 
of the field of the slab on y. This approximation could be relaxed with 
little additional complication. 
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Using the above approximation, we obtain for the field of the H i 
mode (rod) at r = a, from eqs. (55), (56), and (57), 

H. = (iuno^giJoigid), (65a) 

E* = pTHupJB,, (65b) 

P = rMniO-Wlfoio), ( 65c > 

where 

g\ ^ «V« -h\= (3.8/o)», 

Pi = hi — <o 2 € ju « w 2 (e — e )no = u 2 . 

For the slab we have, at r = a, from eqs. (49), (50), and (51), 
setting the arbitrary constant E, a equal to unity, h.„ £zl h, and taking 
into account the exp (p,x) dependence of the field below the slab 

H, z = -2(i^ fJ . )- 1 (T/2d) exp [-p.(Z> + a</> 2 /2)], (67a) 

E.y = -vrHt*n H, z , (67b) 

P. = 2h t {o ) y. )- l dL y , (67c) 

with 

h,tthi£zi kn, 

p.«pi«« = /e(n 2 - 1)*, (68) 

7r/2d = 2.4/o. 

The coupling coefficient C 2 is c 2 /PP„ where 

c = af +T [#*#„ - E. y cos (0)ff,]o>. (69) 

From eqs. (65) and (67), it is apparent that the two terms in the 
integrand in eq. (69) are equal and add up if we make the approxima- 
tion cos & 1. Thus, 

c « 2aE4, f °° H„d<t>. (70) 

J — oo 

Using eq. (67a) for H it , we have 

f +K H.M = -2(ia> Mo )- 1 (T/2d) exp {-p.D)(2v/p.a)i, (71) 

J— 00 

if we make use of the identity 

f +C ° e- bzl dx = (»/&)*. (72) 

Thus, 

c = 4apr 1 0i(i"/i o )- 1 (T/2aVo(0ia) exp (-p.0)(2ir/p.a)* > (73) 
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and 

C 2 = (32/7r) & ?(ir/2ri) 3 (w 3 /i 2 aL 1 ,)- 1 exp (-2uD). (74) 

Since the mode number density is given by eq. (40), the loss 

£ = h(h, z /h 8V )C*L y (75) 

is finally obtained from eqs. (74), (62), (66), and (68), 

£ = 340n- 1 (n 2 - l)-*^ 5 ) -1 exp [-2(n 2 - l)*fcD]. (76) 

The loss in dB/km is obtained by multiplying the r.h.s. of eq. (76) 
by 8.7 X 10 9 , the /xm being used as the unit of length. Thus, for 
n = 1.41 and n = 1.01 we have, respectively, 



£dB/km = 1.35 X 10 9 X^(a/X)- 5 exp (-12.5D/X), 
£dB/km = 675 X 10 9 X^(a/\)- 5 exp (-1.76.D/X), 



n = 1.41, (77) 
n = 1.01. (78) 



For example, if D = 0.15 /xm, n = 1.41, X = 1 nm and a = 40 itm, 
we find that the radiation loss of the Hoi mode is £ = 2 dB/km. If 
D = 1 pm, n = 1.01, X = 1 itm, and a = 40 /im, the loss is as high 
as 1140 dB/km. The radiation loss is shown as a function of a/X and 
D/\ in Figs. 5 and 6 for a wavelength of 1 /mi, and for n = 1.41 and 
1.01, respectively. The amount of loss required to prevent the power 




Fig. 5 — Radiation loss in dB/km of the rod Hoi mode in the slab as a function of 
spacing D with the rod radius a as a parameter, for n r0 d = n 8 ub = 1.41. These curves 
are valid for large values of D. 
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Fig. 6 — Continuation of Fig. 5 f or n = 1.01. 

transferred to the H i mode to be transferred back to the HEu mode 
and to cause pulse spreading depends on the fiber irregularities and is 
not accurately known. 

The above results are approximate and, to some extent, incomplete. 
In particular, the perturbation method that we used is not accurate 
when D is small. Also it would be useful to ascertain that the radia- 
tion losses of the other higher-order modes are at least equal to the 
loss calculated for the Hoi mode. For some of these higher-order modes 
of the rod, it is necessary to take into account the higher-order modes 
of the substrate, both E and H, and this involves some complication. 9 
In spite of these limitations, our result, eq. (76), should provide pre- 
liminary information concerning the mode-selection mechanism 
afforded by 2-dimensional mode sinks. In particular, the very fast 
dependence of the loss on the rod radius (a -5 ) indicates that very 
large rods cannot be used if single-mode operation is to be achieved 
in air. However, if the gap between the rod and the slab is filled up 
with a material whose permittivity is only slightly smaller than the 
rod and slab permittivities, the rod radius a and the spacing D can be 
large, as Fig. 6 suggests. 
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APPENDIX 

Limit forms of the propagation constants in optical fibers 

Two approximations can be made, applicable to low-order modes 
in highly multimoded fibers and to fibers with small transverse varia- 
tion of permittivity. A simplified presentation is given in this appendix. 

Low-order modes propagating in highly multimoded fibers cor- 
respond to waves propagating almost along the axial direction, z. 
The propagation constant h is therefore close to kn if n denotes the 
refractive index on axis. If the fiber refractive index is a constant 
within some contour and assumes a lower value outside that contour, 
the wave near a section of the contour can be assumed plane. Because 
it is incident at grazing angles, the electric and magnetic fields tend to 
zero compared with their values in the bulk. Thus, the electric and 
magnetic fields at the boundary of a dielectric rod vanish, compared 
to their values on axis, as the transverse dimensions of the rod tend 
to infinity for a given mode number. 

For a round fiber with refractive index n and radius a, the exact 
equation defining h is, using the notation of the main text, 6 

T nV;(tti) + KIM ir JIM + Kl(u 2 ) i 



VJU 
J L u iJ> 



UiJ„(Ui) UiK v (Ut) J[ U\J v {Ux) UiKriUz) J 

= [1^1 m 

the axial and azimuthal variations of the field being denoted 
exp (ihz + iv<t>) and 

Mi = ga = (k 2 n 2 — h 2 )m, 

u 2 = pa= (hr - fc 2 )*o, (80) 

k = 0)(e o n )K 

In the limit a — > « , u-> tends to infinity and the second terms in the 
brackets on the l.h.s. of eq. (79) vanish [Kl{x)/K,(x) — * — 1 if x — * *>]. 
On the r.h.s. of eq. (79), h can be replaced by kn. Thus, it is apparent 
that eq. (79) becomes 

Jl(ui)/J„(ui) = ± v/ui, (81) 

or, equivalently, using well-known formulas involving Bessel's func- 
tions and their derivatives : * 

J^ti(wi) = 0. (82) 



* We have vJ, =F xJ' v = xJ, ± i and (for later use) vK u ± xK' v = =FxK, ± i. 
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Table I 



,-2-10 1 2 

EH: J- 3 = J_, = J-i = Jo = Ji = 0: HE~| 

v<0<{ (M--2) (M=-D (m = 0) ( m = 1) (M-2) lr>0 

.HE: J_i = 7 = Ji = J 2 = J 3 = 0: EH. 
[Note: J_ = 0«=>/, = 0] 



For symmetry reasons, modes with opposite values of v have the same 
propagation constants. For v = 2, for instance, the propagation con- 
stants of the two sets of modes are given by the roots J\ and J 3 - For 
v = — 2, they are given by the roots J- 3 and J_i. However, these are 
the same because J- v = {-)"J». Equation (82) was given by Snitzer. 7 
For the HEn {v = 1) and H i (v = 0) modes, the relevant solutions of 
eq. (82) are the first roots of 



Jo(ui) = 0, wio = 2.4 



and 



(83a) 
(83b) 



Ji(wi) = 0, wio = 3.8- • •. 

These are the results used in the main text. 

Because the modes Hoi, Eoi, and HE 2 i have almost the same propaga- 
tion constants (see Table I), the validity of the calculations given in 
the main text can be questioned where the mode Hoi was considered 
independently of the two other modes. It is therefore important to 
evaluate the actual splitting between these three modes. For simplicity, 
we consider only the Hoi and Eoi modes. The expressions giving the 
exact propagation constants of the H i and E i modes are, setting v = 
in eq. (79), 

Ji(ui)/uiJo(ui) = —KiM/uiKoim), 



and 



Setting 



where 



Ji(ui)/uiJo(ui) = -n 2 Ki(u 2 )/u2K (u2), 
Wi = u + 5, K<1, 



J\{u ) =0, u = 3. J 
on the l.h.s. of eqs. (84) and (85) and 

w 2 = k(n 2 — l)*a 



(Ho), (84) 

(E ). (85) 

(86) 

(87) 

(88) 
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on the r.h.s., we obtain for the difference Ah in propagation constant 
between the H i and Eoi modes 

Aha = (3.8/kna) 2 (l - l/n 2 )». (89) 

Except for a numerical factor, this result is the same as for a slab (see 
Section 4.1). If a = 10 /xm, n = 1.41, and X = 1 fxm, the beat wave- 
length 2ir/Ah is, from eq. (89), equal to 5 cm. 

The individuality of the H i mode is preserved and the calculations 
given in the main text are valid if the loss £ is small over that length 
(e.g., £ « 1 dB/cm for a = 10 /xm). In fact, this restriction on £ 
may be even less stringent than that calculated above because the de- 
generacy between the three modes may be lifted further by the pres- 
ence of the slab when the coupling is increased. 

The second approximation referred to at the beginning of this ap- 
pendix is the scalar approximation widely used in optics. If the trans- 
verse variations of the medium permittivity are small, the x and y 
components of the field satisfy approximately the scalar Helmholtz 
equation 

(d*/dx* + d*/dy*)E x + [&n*(x, y) - h 2 ~]E x = 0. (90) 

A similar equation holds for E y , which need not be written down. 

Because all quantities are bounded in eq. (90), E x and its first 
derivatives are continuous functions of x and y. 

For the rod considered earlier, eq. (90) becomes, assuming an 
exp (in4>) dependence of E x on 0, 

d 2 E x /dr 2 + r-'dEJdr + (IbV - h 2 - n 2 /r 2 )E x = 0, r < a, 
d*E x /dr* + r~'dE x /dr + (k 2 - h 2 - n 2 /r 2 )E x = 0, r > a. (91) 

These are differential equations for Bessel functions. The bounded 
solutions of eq. (91) are 



(92) 



E x = JM, g 2 m khi 2 - h 2 , r < a 

E s = AK^ipr), p 2 = h 2 - k 2 , r > a. 

Continuity of E x and dE x /dr imposes 

J^M/K^ut) = (ui/u2)J^i)/KM, (93) 
or, using the transformation formulas given before, 

UiJ h-i(wi) I J m(wi) = u 2 K tl+1 (u 2 )/K ll (u2) ) (94) 

a result previously derived by Snyder 5 from the exact equation, eq. 

TRANSVERSE COUPLING IN FIBER OPTICS II 695 



(79). In the limit a — > «, eq. (94) reduces to 

J,M = 0, (95) 

in agreement with eq. (82). To each value of /z we must associate modes 
corresponding to the two states of polarization of the electromagnetic 
field. This is illustrated in Table I. 

The physical significance of the scalar approximation is that if, 
for instance, a linearly polarized field, solution of eq. (90), is launched 
into a fiber, this field configuration is approximately maintained over 
a certain length. Eventually, however, the polarization is transformed 
because the two electromagnetic modes have slightly different real 
propagation constants as we have seen (for a report of experimental 
observations, see Ref. 8 in which the mode n = ±1 is illustrated in 
Figs. 3 and 4d) and/or different losses. The scalar approximation is 
useful to obtain approximate expressions for the propagation con- 
stants. This approximation is not applicable to the evaluation of 
radiation losses if these losses are polarization dependent. This is the 
case, for instance, if the propagation constant of the rod mode lies 
between the propagation constants of the slab E and H modes. Because 
the split between these two modes is very small, this is unlikely to 
happen unless the optical waveguide has been specially designed for 
that purpose. In that sense, the scalar approximation may be applied 
to problems of radiation losses. 
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